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Abstract 

The out-of-equilibrium dynamics of the 0{N + 1) nonhnear a-model in 1 + 1 
dimensions is investigated in the large N Umit. Regarding the nonUnearity 
as the effect of a suitable large coupling limit of the 0{N + 1) (j)^ model, we 
first of all verify that the two limits commute, so that the 0(oo) nonlinear a 
model is uniquely defined. Such model can be completely renormalized also in 
the out-of-equilibrium context, allowing us to study the consequences of its 
asymptotic freedom on the time evolution far from equilibrium. In particular 
we numerically study the spectrum of produced particles during the relaxation 
of an initial condensate and find no evidence for parametric resonance, a result 
that is consistent with the presence of the nonlinear contraint. Only a weak 
nonlinear resonance at late times is observed. 
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I. INTRODUCTION AND SUMMARY 



Recently, a great deal of attention has been paid to the relaxation of a bosonic condensate 
in interaction with its quantum and/or thermal fluctuations. Some of the main results of this 
research program have been obtained in the description of the inflationary dynamics, where 
one has to consider the expectation value of a scalar field rolling down towards equilibrium. 
Since similar scalar theories can be used to describe the low energy features of hadronic 
physics (like the interaction among pions), this subject has been studied also in connection 
with the formation of Disoriented Chiral Condensates (DCC), which may happen after the 
collision of two ultrarelativistic nuclei. 

Thus, much work has been done about the quantum evolution out of equilibrium of the 
model in 3 + 1 dimensions As is well known [0], the renormalized theory is trivial. 

Practically, this means that we should consider the model as an effective theory, keeping the 
ultraviolet cut-off A much smaller than some Landau scale. The logarithmic dependence on 
A should disappear from the renormalized quantities, while a weak inverse power dependence 
remains. 

If we want to push the application of non equilibrium techniques to more fundamental 
theories, like QCD, we should consider that the ultraviolet properties change drastically. In 
those cases, in fact, there is no Landau Pole in the ultraviolet and the renormalized cou- 
pling becomes smaller and smaller as the momentum scale increases. This corresponds to 
the property of asymptotic freedom, whose presence justifies self-consistently the perturba- 
tive renormalization procedure and allows in principle to perform the infinite cut-off limit 
smoothly. 

Motivated by this consideration and by its intrinsic relevance both in Quantum Field 
Theory and in Statistical Mechanics, we analyse in this paper the dynamical properties of 
the nonlinear a model in 1 + 1 dimensions, which is asymptotically free in the ultraviolet 
@]. Thus, A can be pushed to infinity rigorously and there should exist a renormalized 
out-of-equilibrium dynamics, completely independent of the ultraviolet cut-off. 

The linear and nonlinear a models in 3 + 1 dimensions were introduced in elementary 
particle theory in order to provide a useful model of the low-energy strong interaction sector, 
which was able to realize the SU{2) x SU (2) current algebra and the Partial Conservation of 
Axial Current (PCAC) and satisfy the corresponding low energy theorems 0. Afterwards, 
the nonlinear a model has been considered fruitfully in many areas of Quantum Field Theory 
and Statistical Mechanics, mainly in the description of quantum spin chains and 2D spin 
models [§] and, quite recently, of disordered conductors and of quantum chaos . 

We consider in this paper the 0{N) invariant model, in the large N limit. First of all, we 
present in section y an original derivation of the relevant result that the large N limit and 
the large coupling limit, which turns the linear model in the nonlinear one, commute. More 
precisely, we show that the same classical (in the sense of Yaffe [|TT|) hamiltonian, which 
describes the quantum dynamics of the model in the large approximation, is obtained, 
no matter in which order we perform the two limits. Previous studies on this subject were 



presented in |]T^, using perturbative techniques and the derivative expansion for the model 
in 3-1-1 dimensions, with the conclusion that the divergent terms are universal, while finite 
parts do differ when taking the large coupling limit on the quantum corrections on the linear 
model, or calculating the same quantum corrections on the nonlinear model. In our case we 
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find instead that the large coupling limit of the 0(oo) 0^ model is completely equivalent to 
the large limit of the quantum nonlinear a model. We also derive the evolution equations 
for the nonlinear model at the leading order in the expansion, in the case of an initial 
field condensate different from zero. We implement the constraint by the use of a Lagrange 
multiplier, which we denote m^, since it enters the dynamics as a squared mass. We show 
that the usual renormalization procedure, which makes the bare coupling constant depend 
on the UV cutoff, is sufficient to get properly renormalized, that is UV finite, evolution 
equations. Moreover, we characterize the ground state of the model, giving an interpretation 
of the dynamical generation of mass (the so-called dimensional transmutation) in terms of a 
compromise between energetic requirements and the constraint. We conclude by describing 
suitable initial conditions for the condensate and the quantum fluctuations. We want to 
emphasize that, while the approach we follow in the study of the dynamical evolution in 
this field theory has become by now quite standard and is very similar to that of ref.s 
the different dynamical properties of the nonlinear a model leads to different and new 
results, most notably the absence of parametric resonances and the consequent need for a 
detailed study of the scaling properties with respect to variation of A, which we describe in 
the second part of this work. 

In sec. |ITT| we present the analysis of the numerical evolution for the condensate and the 
Lagrange multiplier as well as for the number of particles created during the relaxation of 
the condensate (the quantum fluctuations). Remarkably, we do not find any period of expo- 
nential growth for the fluctuations. Actually, no spinodal instabihties were to be expected, 
since the symmetry is always unbroken in 1+1 dimension. But there occurs also no para- 
metric resonance, as takes place instead in the unbroken symmetry scenario of the large— A^ 
0^ model in 3 + 1 dimension. This is due to the quite different nonlinearities of the a— model 
and in particular to the nonlinear constraint [see eq. (2.29)] which sets an upper bound 
to the quantum infrared fluctuations [see fig. In fact, even if the constraint disappears 
as the bare coupling constant Ab vanishes in the infinite UV cutoff limit (asymptotic free- 
dom), the quantum fluctuations in any given finite range of momentum remain constrained 
to finite values, as implied by the possibility of fully renormalize the model, including the 
constraint [see eq. ( |2.46| )]. Because of this and of the reduced momentum phase space, we 
observe that the damping of the condensate is not as efficient as in the large— A^ 0^ model 
in 3 + 1 dimension with unbroken symmetry. As a matter of fact our data do not even allow 
to establish for sure that the condensate will eventually relax to zero (in the case with zero 
angular momentum, see below). 

In particular, we carefully study, by numerical fit and self-consistent analytic computa- 
tions, the asymptotic evolution of the Lagrange multiplier. The estimated dependence of 
its asymptotic value, m{ooY, on the initial condensate po, turns out to be very well approx- 
imated by an exponential, which is the exact dependence of m(0)^ [at infinite UV cutoff, 
see eqs. ( |2.36| ) and ( |2.37| )]; remarkably however, the prefactor in the exponent is changed 
[see eq. (|3.?1 )]. We also verify that, after the proper renormalization, the dependence of the 
asymptotic values on the UV cutoff A is only by inverse powers. As far as the emission of 
particles is concerned, we considered three different reference states: the initial state, the 
adiabatic vacuum state and the equilibrium vacuum state, that is the true ground state of 
the theory. The numerical results suggest a weak nonlinear resonance, yielding a relaxation 
of the condensate via particle production driven by power laws with non universal anoma- 
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lous exponents, a result similar to what found in 0] for the asymptotic dynamics of 0^ in 
3 + 1 dimensions. However, more numerical as well as analytical work is necessary for a 
better quantitative estimates. 

Finally, since we allow the condensate to have a number n of components larger than 1, 
we are able to study the evolution of configurations with non-zero angular momentum £ in 
the internal space of the field [see eq. (2.27)]. In this case we find numerical evidence for 
an adiabatic spectrum broader than in the case £ = [see figure ITT]], suggesting a stronger 
coupling with hard modes. Again, more work is necessary for a better understanding of this 
issue. 

There are, of course, physically relevant issues which are not addressed in this paper, 
like the effects of sub leading (in terms, which may lead to genuine thermalization, and 
the consequent non-uniformities between large time and large N , or the comparison with 
the evolution of spatial inhomogeneous condensate [^ . 

It nevertheless appears evident that, when compared to the unconstrained 0^ model, the 
dynamics of this constrained model has qualitatively different properties, which deserve a 
separated and more detailed analysis. 



II. THE 0(oo) NONLINEAR a MODEL IN 1 + 1 DIMENSIONS 



A. Definitions 



The classical Lagrangian of the 0{N + 1) a model is given by 



(2.1) 



where is a multiplet transforming under the fundamental representation of 0{N + 1) and 
constrained to the A^— dimensional sphere of radius A~^/^: 



0-(/) = 1/A 



(2.2) 



A may be regarded as the coupling constant, since the sphere flattens out in the A — >^ limit. 
The Hamiltonian corresponding to (|2.1|) reads 



H = - dx 
2 



J' + (40)' 



i<j 



(2.3) 



where Jtj = (piUj — 0j7rj is the angular momentum on the sphere, tTj being the momentum 
conjugated to 0^. This Hamiltonian can also be obtained as the g oo limit of the linear 
model 



dx 



(2.4) 



where (f) is now unconstrained and the potential V may be taken of the form 



v{u) = ^-{u-i/xr 



(2.5) 
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The quantum version of the Unear model defines a textbook Quantum Field Theory (apart 
from the nontrivial strong coupling limit g oo). The quantum version of the nonlinear 
model ( |2.3| ) may be written instead 



(2.6) 



1 + Abtt^ 



H =-J dx [-A+oo^ (d^a)^ 
where we have used the projective coordinates (ai, ... , oat) on the sphere, namely 

1/2 

4>j = UJOj , Ab (j)N+l = UJ -1 , UJ 

SO that 

and the 0{N + 1)— symmetric functional Laplacian reads 

Aix)=ujix)-''—^ujixf~^ ^ 



6aj{x) 



6a j (x) 



(2.7) 



(2.8) 



(2.9) 



while a in eq. ( p.6|) is a multiplicative operator. We have replaced the coupling constant A 
with Ab (the bare coupling constant) to stress the fact that in Quantum Field Theory it is 
generally cut-off dependent. 



B. The iV ^ oo limit 



Now we derive the quantum dynamics in the large N limit, applying a general technique 
already used in the analysis of the 0^ dynamics in finite volume |]13| and based on well- 
known work by Yaffe |Tl|. If we consider the nonlinear model as a limit of the 0^ linear 



model (being this true at least at the classical level), we have to take two limits and we might 
wonder whether it is legitimate to interchange their order. To be more specific, if we first 
perform the large N limit in the linear model, we get a classical dependent unconstrained 
Hamiltonian that admits a definite nonlinear limit H°° as (? ^ cxo. We verify here that 
indeed the same Hamiltonian follows if we start directly form the nonlinear quantum 
Hamiltonian ( |2.6|) and take the N oo kla. Yaffe. 

Consider the quantum Hamiltonian of the linear model, with the couplings suitably 
rescaled to allow the large N limit 



ifi, = dx 



V{u) = ^{u-N/Kf (2.10) 



According to (a slight extension of) Yaffe's rules, the quantum dynamics described by the 
N ^ oo limit of the model is described by a classical Hamiltonian, which is the large N limit 
of the expectation value of the quantum hamiltonian ( p. 10 ) on a set of generalize coherent 



states, labelled by the parameters defined in eq. ( p.l2[ ). We end up with the following 
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classical hamiltonian 



dx 



+ - y dx dx" v{x, x')w{x', x")v{x", x) 



^-n^ + ^{d.cP)' + V{cP' + w{x,x)) 



+ / dx 



^{x,x) - ^dlw{x,x'] 



V{u) 



(2.11) 



Tr{x) 



lim — 



(2.12) 



(2.13) 



where the classical canonical variables are defined as 

/ v^(0(x)) 
yiV(7r(x)) 

W{X, X'} I 7V^<x) 7V (0(x) ■ 0(x'))con: 

and the nonvanishing Poisson brackets read 

{(t)j{x) , vrfe(x')}p B = SjJix - x') 
{w{x, y) , v{x', yOlp.B. = Six- x')6{y - y') + 6{x - y')5{y - x') 

It is understood that the dimensionality of the vectors 4>{x) and 7r(x) is arbitrary but finite 
[that is, only a finite number, say n + 1, of pairs {(f>{x) ,Tt{x)) may take a nonvanishing 
expectation value as — oo]. Thus, the index j may run form 1 to n + 1, where n + 1 is 
the number of field components with non zero expectation value. 

The g ^ oo limit on the classical Hamiltonian is straightforward and reintroduces 
the spherical constraint in the new form 

n+l 

^ ^ diag(w;) = 1/Ab (2.14) 
whose conservation in time implies 

(f) ■ 77 + dia.g{wv) = (2.15) 
where we have introduced the condensed notation 



ab){x,y) = / dz a{x, z)h{z,y) , diag(a)(x) = a(x,x) 



(2.16) 

Let us now come back to the quantum Hamiltonian ( |2.6| ) of the non-linear model. First of 
all we perform a similitude transformation of the Laplacian, to cast it in a form suitable for 
the application of Yaffe's method: 



- A 



1 2\a^ 

2 ^ ivj iV 



N 



+ 1 



(2.17) 



where 6ij{x) is the obvious multiphcation operator and i3j{x) = —i6/6aj{x) its conjugated 
momentum. Now, after the rescaling Ab X^/N in eqs. (p.7|), by the usual rules in the 
N oo limit we obtain the classical Hamiltonian 



1 

"2 



+ diag(xr/x) + ^diag(r7 ^) 



dx 

(9^.a)^ + dA'il{x,x') 



(2.18) 
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where 



n 



1 + Ab[a2 + diag{r])] 



(2.19) 



and, just as in eq. ( p.l2D , 



f3{x) 
ri{x, x' 
\x{x,x')J 



lim — 



/ VN{a{x)) 
VN{hx)) 
(q;(x) ■ 6l{x'))c 
V(^(a;)-^(x'))o 



(2.20) 



are classical canonically conjugated pairs, with Poisson brackets identical to those in eq. 
( p.l3| ). We take the indices of the classical fields cx and (3 to run form 1 to n, having assumed 
that only the first n components of their quantum counterparts may have expectation values 
of order \fN . 



To show that the classical Hamiltonian is equivalent to the g —>■ oo limit of H^, we 
need only to solve the spherical constraint ( 2.14 ) that emerges in that limit. This amounts 
to the canonical parameterization of the constrained pairs {4>, tt) and {w, v) in terms of the 
projective ones (a, /3) and {ri,x)- It reads 



\i/2n 



(Pn+1 



w{x,x') = Q{x)Q{x') ri{x,x') , v{x,x') 



n-1 

-a - f3~ diag(?7x) 
xix,x') 6{x — x') 



n{x)n{x'] 

and in particular it implies, besides ( p. 14]) and (|2.15|) , 



+ 



n{x) 



(2.21) 
(2.22) 

7r„+i(x) (2.23) 



TT^ + diag{vwv) 



{d^(t)f + dA'W^x.x') 



P"^ + diag(x??x) 
(9^.a)^ + d^d^:r]{x,x'] 



(2.24) 



This result proves the complete equivalence between the g oo limit on the leading 1/N 
term of the linear model (which imposes the new spherical constraint) and the N ^ oo limit 
of the quantum model directly formulated on the constraint manifold. 

Before closing this section, it should be noticed that, even though we gave the basic 
definitions and performed the entire computation for a field theory in 1 + 1 dimensions, the 
results in sections [11 A| and [11 B| remain valid also for a. {D + 1)— dimensional theory, the only 
change being in the dimensionality of the integrals. 



C. Dynamical Evolution 

Let us now derive the evolution equation for this system in the case the field has a non 
zero, albeit uniform, expectation value 4> in the initial state. The 2— point functions depend 
only on the difference x — x', and can be parametrized by time-dependent widths 0"^: 

wix, x') = ^ al e"'^^-^'^ , vix, x') = ^K--^') (2.25) 
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where A is the ultraviolet cut-off. 

In this case of translation invariance, in practice one can always take n = 1 owing to the 
0{n + 1) symmetry of H^. Thus, we choose to have only two non-zero components. In 
other words the condensate will move on the plane specified by the initial conditions for (f) 
and its velocity. Using eq. ( p.25| ), we may write the Lagrangian density corresponding to 
the {g oo limit) of the Hamiltonian (|2.11[ 



as 




1 



(2.26) 



We have kept into account the constraint by introducing the Lagrange multiplier m^. The 
corresponding Euler-Lagrange evolution equations read, in polar coordinates 

p + m^p- — = (2.27) 
at + {k'' + m^) dfc - ^ = (2.28) 
p2 + S - ^ = (2.29) 

with the definitions £ = p^9 (the conserved angular momentum of the condensate) and 

S = diag(.^) = i^ol (2.30) 

The first thing we can do is to look for the minimum of the Hamiltonian, that is the ground 
state of the theory which corresponds to the vanishing p, p, cta;, at and £. The equations to 
solve are: 

rn?p = 

2 , 2^ _ 1 . (2-31) 



{k +m )ak 



^4 



The solution m = is not acceptable, because it yields a massless spectrum for the fluc- 
tuations and gives an infrared divergence that violates the constraint. This is nothing else 
than a different formulation of the well-known Mermin-Wagner-Coleman theorem stating 
the impossibility of the spontaneous symmetry breaking in 1 + 1 dimensions . Thus, the 
unique solution is: p = and = ^{k'^ + m^)~^/^. 

This result allows for an interpretation of the mechanism of dynamical generation of 
mass as the competition between the energy and the constraint: in order to minimize the 
"Heisenberg" term in the Hamiltonian, the zero mode width, that is do, should be as large 
as possible; on the other hand, it cannot be greater than a certain value, because it must 
also satisfy the constraint. The compromise generates a mass term, the same for all modes, 
which we call rrieq 
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We can take the mass at equilibrium as an independent mass scale defining the theory, 
as dictated by the dimensional transmutation, and the relation between this mass scale and 
the bare coupling constant is read directly from the constraint (2.29) 



Ab 27r 



log 





1 \ 











1 , 2A ^ 
— log + 

ZTT rrir. 



m. 



cq 



"cq 



A2 



(2.32) 



When the system is out of equilibrium, the Lagrange multiplier m may depend on time. 
Its behavior is determined by the fact that the dynamical variables must satisfy the con- 
straint. After some algebra, this parameter can be written as: 



m 



Ab ( + ^ + e 



e 



A dk 

-A 2^ 



(2.33) 



We can describe the quantum fiuctuations also by complex mode functions z^, which are 
related to the real function ak by: 



4 Ok 



1 

2 



\Zk\ 



(2.34) 



One can recognize in the second term on the r.h.s. of the last equation in ( p. 34 ) the cen- 
trifugal energy induced by Heisenberg uncertainty principle. 

We choose the following initial conditions for this complex mode functions: 



Zk{0) = -^/== , 4(0) = -i^JuJk/2 



(2.35) 



where ook = \/k'^ ^ o? and a is an initial mass scale. It is worth noticing here that such a 
form for the initial spectrum of the quantum fiuctuations does not allow for an initial radial 
speed for the condensate degrees of freedom, unless we start from po = 0. This is easily seen 
by differentiating (2.29) with respect to time. 

Moreover, we should stress that a might be different from the initial value of the Lagrange 
multiplier. In fact, once the initial value for p is fixed, a can be determined by means of the 
constraint equation and it turns out to be 



a(po 



cq 



eq 



exp (27rp2) |i 
exp (27rpo) 



2 

1 + 



1 + 11 + 



eq 



A2 



+ exp(47rpo) 



'1 + 



eq 



A2 



eq 



A2 



On the other hand, the initial value for the Lagrange multiplier is given by 

/.2 /2 \ 

2 1 \ \ 22i 

mo = a + Ab I Po +^ - " Po I 



(2.36) 



(2.37) 



that is equal to the initial mass scale o? only if we push the ultraviolet cut-off to infinity. 
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To properly control for any time the ultraviolet behavior of the integrals in eqs. (p.30|) 
and (p. 33 ), one should perform a WKB analysis ||l^ of the solution. One finds the following 
asymptotics for the mode functions: 



2k Jo 



Zk{t) = ZkiO) exp \ikt-— / dt' m\t') - — m\t) - m^(0) 



(2.38) 



From the above formula it is clear that the logarithmic ultraviolet divergence in E is com- 
pletely determined by the initial spectrum. For the divergent integral in eq. ( |2.33| ) the 
situation is more involved. Explicitly one finds: 



J-A Zn Zn /X 



and 



where 



A dk 



— (|4(t)|2 - e\zk{t)\') = mitf S(0) + GfW 



A 27r 



(2.39) 



(2.40) 



SF(/i;t) 



A dk 

-A 2^ 

A dk 



\Zk{t)\ 



2 o{\k\-^^) 



2\k\ 



(2.41) 



0F(t)=/ TT \zk{t)\' - k'\zk{t)\' - m{tY\zkm 

K Z-K ^ 



have finite limits as A — cxo. We have introduced in the above formulae a subtraction point 
[X. There correspond a renormalized coupling constant A running with /i, as the A — oo 
limit of the relation 



1 , A 1 1 

log ^^TT + 



27r "/i Ab(A) A(/i) 







and a renormalized constraint 



p(t)2 + SF(/i;t) 



1 



A(/i) 







(2.42) 



(2.43) 



With this definitions, the equilibrium mass scale can be written as 

27r 



m. 



eq 



2/iexp 



which by consistency with eq. ( 2.36|) implies 

A(a/2) = \IpI 

In conclusion we can rewrite the constraint and the Lagrange multiplier as 



+ Sf(/u) - 



-— = , m = — 
Po 



p2 + — + Gf 
p2 



(2.44) 



(2.45) 



(2.46) 
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For large but finite UV cutoff tfiese expressions retain a inverse power corrections in A. In 
tlie actual numerical calculations whose results will be presented in the following section, 
we used the "bare" counterparts of eqs ( |2.46| ) with finite cutoff and the definition ( |2.32| ) of 
the bare coupling constant is used to reduce to inverse power the cutoff dependence. 

Let us conclude this section by summarizing the steps we need to do, before trying to 
solve numerically the equations of motion. Once we have fixed the UV cutoff, the equilibrium 
mass scale rrieq and the initial value for the condensate po, we can determine the initial mass 
scale a in the fluctuation spectrum from eq. ( p.36| ), which in turn gives the initial conditions 
for the complex mode functions [cfr. eq. ( |2.35 )]. Now, we need to specify the remaining 
initial values for the condensate, namely its velocity po and its angular momentum i, which 
must be consistent with the constraint ( p.l5|) . Finally, eq. ( p.46| ) completely determines the 
initial value for the Lagrange multiplier mo, which has exactly the same infinite cutoff limit 
as a, but differs significantly from it for finite cutoffs. 



III. NUMERICAL RESULTS 



We have studied numerically the following evolution equations 



+ m 







m 



m Zk 



Ab 



^ dk 



(3.1) 



UK / 

9^ 



A 2n 



where = 0i + 202 = pe^^, 9 = i and |0p = + £^/p^, while the bare coupling constant 
Ab is given by eq. ( p.32| ). The initial conditions for 0, and Zk [see eq.s ( p.35| )] must satisfy 
the constraints ( 2.14 ) and (|2.15 ), that are then preserved by dynamics. 

In the classical limit the quantum fluctuations Zk disappear from the dynamics. In that 
case the stationary solutions are trivial: 



Pit) = A, 



-1/2 



(3.2) 



with arbitrary value for the angular momentum i. Thus there are stationary solutions 
corresponding to circular motion with constant angular velocity. 

When we include the coupling with quantum fluctuations, we still obtain stationary so- 
lutions, parametrized by i which assumes arbitrary positive values. They have the following 
form: 



Pit) 



\ 



m{t) = nieqX 



(3.3) 



where x depends on i through 

2ni 



+ sinh ^ 



A 



Tn^qX 



sinh ^ 



A 



m. 



(3.4) 



which reduces to xlogx = 27r£/meq in the inflnite cut-off limit. 
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A. Evolution of condensate and Lagrange multiplier 

In order to control the dependence of the dynamics on the uhraviolet cutoff, we solved the 
equations of motion for values of A ranging from Smgq to 20meq, with an initial condensate 
ranging from po = 0.2 to po = 0.7. We mainly considered the case ^ = 0. A typical example 
of the time evolution of the relevant variables is showed in Fig.s ^ and |^. Figure ^ shows 
the evolution of the Lagrange multiplier m(t)^ for A/meq = 20; in this case, its starting 
value is 2.630632. Due to the lack of massless particles, the damping of the oscillations of p 
and is very slow, as already noticed in 0] for the linear model in 1 + 1 dimensions; the 
dissipation is not as efficient as for the unbroken symmetry scenario in 3 + 1 dimensions, 
because of the reduced phase space. A detailed numerical study of the asymptotic behavior 
and a FFT analysis of the evolution allows a precise determination of the asymptotic value 
and the main frequency of oscillation of the Lagrange multiplier: 

m{tr=ml + Pf + 0[^) (3.5) 
where the function p(t) turns out to be 

~ Acos(2moot + 7ilog'': + 72) (3.6) 

The logarithmic dependence in the phase could be justified by self-consistent requirements 
(see below), along the same lines of the detailed calculations performed in ref. 0] in a similar 
context. Numerically it is very difficult to extract and we do not attempt it here. Comparing 
further our result with that reported in ref. we should emphasize that we do not find 
any oscillatory component of frequency 2mo, as happens instead for the 0^ model in 3 + 1 
dimensions. Moreover, as figure § shows, both the asymptotic mass moo and the amplitude 
A depend on the ultraviolet cutoff A. This dependence may be fitted with great accuracy 
through a low order polynomial in 1/A^, showing that the standard renormalization holds 
at any time, as anticipated by the WKB analysis. Therefore, the extrapolated parameters 
and A give us information on the fully renormalized physical theory (in the large N 
approximation). The table below collects the values of for different values of A and 
of the initial condensate po. The last column contain the extrapolation to infinite cutoff, 
obtained by the low order polynomial fit. The empty cells in the last row correspond to a 
UV cutoff so small that the exact turns out to be negative; these values are excluded 
from the fit. 



Po 


A = 5 


A = 6 


A = 7 


A = 8 


A = 9 


A = 10 


A = 11 


A = 12 


A = 13 


0.2 


1.3073 


1.3047 


1.3032 


1.3022 


1.3014 


1.3010 


1.3006 


1.3004 


1.3001 


0.3 


1.8888 


1.8766 


1.8693 


1.8646 


1.8614 


1.8591 


1.8574 


1.8561 


1.8551 


0.4 


3.3869 


3.3162 


3.2747 


3.2482 


3.2303 


3.2175 


3.2082 


3.2011 


3.1956 


0.5 


8.7094 


7.9915 


7.6082 


7.3764 


7.2246 


7.1193 


7.0432 


6.9861 


6.9424 


0.6 


206.03 


52.433 


35.564 


29.2276 


25.969 


24.016 


22.732 


21.835 


21.178 


0.7 


















238.12 
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Po 


A = 14 


A = 15 


A = 16 


A = 17 


A = 18 


A = 19 


A = 20 


A = oo 


0.2 


1.300 


1.2998 


1.2997 


1.2996 


1.2995 


1.2995 


1.2994 


1.2989 


0.3 


1.8543 


1.8536 


1.8531 


1.8527 


1.8523 


1.8520 


1.8517 


1.8493 


0.4 


3.1912 


3.1877 


3.1848 


3.1824 


3.1805 


3.1788 


3.1774 


3.1643 


0.5 


6.9080 


6.8804 


6.8580 


6.8395 


6.8241 


6.8111 


6.8001 


6.6990 


0.6 


20.682 


20.295 


19.988 


19.740 


19.536 


19.366 


19.223 


16.964 


0.7 


177.645 


146.523 


127.66 


115.07 


106.11 


99.442 


94.302 


68.207 



A similar table can be provided for the amplitude A in the eq. ( p.6|) . The values 
extrapolated to infinite cutoff in a similar fashion as before, turn out to be: 



Po 


0.2 


0.3 


0.4 


0.5 


0.6 


0.7 


A{A = oo) 


0.539 


0.701 


0.924 


1.39 


2.34 


4.30 



However, this fit is not as accurate as that for moo- 

It is interesting to observe that at large UV cutoff has an exponential dependence on 
Po analogous to that of mo (which coincides to a at A = oo). Most remarkably the prefactor 
in the exponent is modified by the time evolution: we find 



mz. ~ exp(27 pi) , 3.5 < 7 < 4.5 (3.7) 



4 



The determination of 7 is rather rough due to the uncertainties in the values of moo extrapo- 
lated to A = 00 at larger po- Notice in any case that the analog of 7 for mo is 27r = 6.28 . . . . 

We also performed some computations for i > 0, with the following results: if we start 
from an out of equilibrium value for p, it will relax through emission of particles towards 
a fixed point, different from the equilibrium value determined by eq. ( |3.3| ). Figures || and 
^ show such a situation for i = 1.0, p(0) = 0.3 and A = lOmeq. In that case we have 
X = 1.000057, while the mean values of the asymptotic oscillations are poo = 0.4203 and 
m^ = 32.0294. 

Before closing this section, we should comment a little further on the evolution of the 



condensate p. When i = 0, fig. shows that the oscillations are actually around zero. 



However, from the available data, it is not possible to decide whether the amplitude will 



eventually vanishes or will tend to a limiting cycle (see fig. p!4D . 

On the other hand, in the case of £ 7^ 0, it is already clear that the condensate does not 
relax to the state of minimum energy compatible with the given value of £, which would 
correspond to the circular orbit with radius given by eq. (|3.3|) . However, it may still relax 
to a circular orbit with a different radius and a different (larger) energy. More detailed 
and longer numerical computations are needed to decide whether the damping reduces the 
oscillation amplitude to zero or not. 



B. Emission spectrum 



Once the evolution equations for the complex mode functions has been solved, it is 
possible to compute the spectrum of the produced particles. First, we should say that the 
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notion of particle number is ambiguous in a time dependent situation. Nevertheless, we 
may give a suitable definition with respect to some particular pointer state. We choose here 
two particular definitions, the same already used in the study of the model 0, plus a 
third one. The first choice corresponds to defining particles with respect to the initial Fock 
vacuum state, the second with respect to the instantaneous adiabatic vacuum state, and 
the third to the equilibrium vacuum (the true vacuum of the theory). The corresponding 
expressions in terms of the complex mode functions are: 



UJk\Zk{t)\^ + 



ad 



Ni 



it) 



UJk 

\zkit)\ 



1 

2 



ad 



UJ 



cq 



1 

2 

1 

2 



ad 



y'P + m(t)2 



(3.8) 



We report our numerical findings on these quantities in figs. |^ - |T^. Since the Lagrange 
multiplier tends asymptotically to a constant value the condensate p{t) oscillates with 
frequency moo and the mode functions Zq{t) with frequency uj{q) = \J + This im- 
plies that particle spectra N^(t) and N^'^i^t) are more and more strongly modulated as time 
elapses, as figs. ^ and |1^ show; on the contrary, N^'^{t) is a slowly varying function of the 
momentum k (cfr. figs. 0-^1), because the oscillations of the mode functions are counter- 
balanced by the time dependence of the adiabatic frequencies \Jk'^ + m(ty. Finally, fig. p!2 
allows for a comparison of the spectra related to different initial values of the condensate. 

Looking at the momentum distribution of the created particles at different times, we 
see the formation of a growing peak corresponding to soft modes. We can give an analytic, 
self-consistent description of this behavior at large times through a perturbative approach, 
similar to the one used in ref. 0. We split the time-dependent Lagrange multiplier in two 
parts, as in equation ( |3.5| ) and we treat the "potential" p{t) / 1 perturbatively, as is done in 
0. We find the following solution: 

sin ujg{t' -t)p{t') 



Age 



iuint 



-ilOnt 



-Zq{t')dt' 



which is equivalent, up to terms of order 0(l/t^), to 



Zg{t) = Ag 



1 + 



Asin^(t) A 
y4sin\I^(t) 



■iUJgt 



Airrioot 



A 

iUgt 



g -\- TTIqq 



+ 



^ujg - mc 



Ug + 



-iiOnt 



Of 4] (3.9) 



with ^{t) = 2moot + 7i logt + 72. The logarithmic dependence is due to the "Coulomb form" 
of the perturbative term p{t)/t in the equations of motion. The expression ( p.9|) displays 



resonant denominators for uj„ 



rrir 



that is g = 0. The perturbative approach is valid as 



long as the first order correction is small compared to zeroth order. Such a condition is 
satisfied if 

A 



At Ug{ujg 



< 1 



rrir 



(3.10) 
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that implies > A/ At for non relativistic modes. Thus the position of the peak found before 
may be interpreted as the result of a weak nonlinear resonance. The asymptotic behavior of 
the condensate and the mode functions related to soft momenta must be obtained through 
non-perturbative techniques, implementing a multitime scales analysis and a dynamical 
resummation of sub- leading terms. A self-consistent justification of the numerical result 
( p.5|) , along with the power law relaxation behavior for the expectation value (with non- 
universal dynamical anomalous dimensions), are likely to be obtained following the line of 
the analysis performed in for the (p^ model in 3 -|- 1 dimensions. 

From the numerical study of the complete spectrum history, we conclude that no expo- 
nentially growing (parametric or spinodal) instabilities are present in the case at hand, as 
apparent from Fig.s ^ - |, which show the spectrum of produced particles with respect to 
the adiabatic vacuum state. 

IV. OUTLOOK 

The natural continuation of this paper is the detailed numerical study of the evolution, in 
order to give a precise picture of the process of mean field dissipation via particle production 
in the framework of this constrained, asymptotically free model. It should be possible to 
determine precisely the power laws that characterize the asymptotic evolution of relevant 
variables, like the condensate, the Lagrange multiplier and the number of created particles. 
After this, one should be able to decide whether, at zero angular momentum, the damping 
leads to the complete dissipation of the energy stored in the condensate or the system evolves 
towards a limit cycle with an asymptotic amplitude different form 0. Also a comparison with 
the linear model in 1 + 1 dimension might be useful to understand the peculiarities of the 
dynamics in a constrained model. 

Moreover, it would be very interesting to study the dependence of the evolution on the 
value of i, the angular momentum of the field in the internal space. As the preliminary 
results presented in this paper show (see figure [Tl|), the asymptotic state is far from the 
state of minimum energy compatible with the given value of £. Remarkably, the adiabatic 
spectrum of produced particles in case of £ 7^ is broader than that one corresponding to 
i = 0, suggesting a stronger couphng with hard modes. 

Most important, the effects of the inclusion of 0{1/N) corrections and the evolution of 
non unigorm condensates should be analyzed also in the framework of the nonlinear a model 
in 1 + 1 dimension. 
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FIG. 3. Evolution of the Lagrange multiplier m'^{t) for A/meq = 20, i = and po = 0.3. In the 
smaller figure there is zoom of the early times. 
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FIG. 4. Asymptotic evolution of Tv?{t) for three different values of the ultraviolet cut-off: from 
top to bottom, A/meq = 5, 10 and 20, £ = and po = 0.3 
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FIG. 5. Evolution of the mean value p for A/meq = 10, po = 0.3 and £ = 1. 
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FIG. 6. Evolution of the squared mass for A/meq = 10, po = 0.3 and £ = 1. 
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FIG. 7. Adiabatic spectrum for trrieq = (solid line), 39.723 (dotted line), 199.006 (dashdot 
line) and 398.11 (dashed line), for A/meq = 20, £ = and po = 0.2. 
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FIG. 8. Adiabatic spectrum for trricq = (solid line), 39.723 (dotted line), 199.006 (dashdot 
line) and 398.11 (dashed line), for A/rrieq = 20, £ = and po = 0.6. 
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FIG. 9. Spectrum with respect to the initial vacuum, for trueq = 39.723 (solid line), 199.006 
(dotted line) and 398.11 (dashdot line), for A/meq = 20, £ = and po = 0.5. 
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FIG. 11. Adiabatic spectrum for trrieq = 0.0 (solid line), 59.895 (dotted line), 119.835 (dashdot 
line) and 179.775 (dashed line), for A/meq = 10, £ = 1 and po = 0.3. 
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FIG. 13. The average value of the condensate p, defined as p = j'^ p{t)dt/T, plotted vs. T, for 
A/mcq = 20, £ = and po = 0.2. 
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FIG. 14. The mean squared fluctuations of the condensate p, defined as J^{p{t) — pY'dt/T , 
plotted vs. 1/T, for the same values of the parameters as in figure [I^. 
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